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Abstract. We investigate the constraints on quintessence arising from both renor- 
malisable and non-renormalisable couplings where the 5d Planck mass is around the 
TeV scale. The quintessence field vacuum expectation value is typically of order the 4d 
Planck mass, while non-renormalisable operators are expected to be suppressed by the 
5d Planck mass. Non-renormalisable operators are therefore important in computing 
the 4d effective quintessence potential. We then study the quantum corrections to the 
quintessence potential due to fermion and graviton loops. The tower of Kaluza-Klcin 
modes competes with the TeV-scale cut-off, altering the graviton contribution to the 
vacuum polarization of quintessence. Nevertheless we show that, as in four dimensions, 
the classical potential is stable to such radiative corrections. 
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1. Introduction 

The case for the existence of a A-like component dominating the currently observable 
universe is now compelling |33] . In the simplest models such a component has two quite 
independent properties. On the one hand, it does not cluster on scales much smaller than 
the Hubble scale, H~ l , and on the other it influences the background evolution of the 
cosmos, causing acceleration at very recent redshifts and giving rise to the coincidence 
problem - why do we appear to be living at a special time in the Universe's history? 
While it is possible to construct models which exhibit only one of these properties, 
eg. (El El, both now have observational support. In particular, the recent detection of 
cross-correlations between the Wilkinson Microwave Anisotropy Probe (WMAP) cosmic 
microwave background (CMB) anisotropics and various tracers of large scale structure 
|Hl EB EH HI] , which are consistent with the decay of perturbations on large scales due 
to an accelerating background, make construction of convincing non-accelerating models 
difficult. 

However, despite its success at the purely phenomenological level, the standard 
ACDM model has almost no deep understanding to back it up. We are in the age of 
precision book-keeping in cosmology, but despite many attempts we do not yet have even 
a well-founded theoretical order-of-magnitude estimate of the size of the cosmological 
constant: most naive field theory calculations disagree by O(10 120 ) with observations, 
yielding perhaps the worst estimate in the history of physics. One can improve the 
situation somewhat by invoking supersymmetry, but it proves generically quite hard to 
construct supergravity vacua with positive cosmological constant 03] . The string theory 
case is even harder [26J. In the absence of any theoretical control over A itself, there is a 
strong temptation to explain the observations by invoking some other mechanism. Some 
proposals utilise the large number of possible string theory vacua, either by appealing 
to the anthropic principle (22j or other quantum effects (2Zj, but a more moderate 
approach is simply to include, among the matter inventory of the universe, some tensile 
matter with appropriate equation of state whose behaviour is under good control. This 
allows us to set A = by supposing that one or more of the string theory proposals for 
cancelling A applies, and then to exclude the complexities of A itself and deal instead 
with the relatively well-understood properties of matter. We will argue that, at least in 
TeV-scale models of quintessence, control is not manifest even in this case. 

Quintessence consists of a scalar field Q, which drives a late-time accelerated 
cosmological expansion via its vacuum expectation value in a rather similar way to 
scalar-field driven inflation. If Q is still rolling today then it must be very light in 
order to satisfy the standard slow-roll conditions and hence its Compton wavelength, 
A c ~ Vqq 2 i IS ver y large (we denote first, second, . . . , Q-derivatives of V as Vq, Vqq, 
. . . , etc.) As a result it only clusters on very large scales, typically greater than 100 
Mpc. 

This is not necessary. In models involving the Albrecht-Skordis potentials, where 
the dark energy reaches a minimum of the potential at non-zero energy, the mass and 
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expectation value of the quintessence is arbitrary, and so the dark energy may cluster on 
all scales after reaching the minimum. Such models are attractive for another important 
reason, for if the quintessence is very light (it typically has a mass tuq ~ 10~ 33 eV), then 
we must find a way to protect this mass from radiative corrections which will otherwise 
spoil the flatness of the potential (see also, eg., Refs. [H1IHI|). This has been studied 
at the one-loop level [T^]. The result depends on which particle species one includes in 
the loops. One typically finds that couplings to bosons are benign [T5] . but couplings 
to fermions are severely constrained. The bounds found by the authors of Ref. are 
extremely stringent and give rise to concern that gravitational couplings alone might be 
strong enough to violate them. Estimates presented in Ref. fH] show that quintessence 
is safe, but this safety is model-dependent and must be assessed carefully. 

In addition, Q must be extremely weakly coupled to standard model fields, 
otherwise it is difficult to see how it could have evaded detection via particle physics 
or cosmological interactions. Despite their overall weakness, such couplings can alter 
standard cosmology in an interesting way [131 El EH], but obtaining them requires 
significant fine-tuning of the renormalisable couplings. 

A more worrying problem is provided by constraints from 4d non-renormalisable 
couplings between the standard model and the quintessence field. Such couplings 
are generically expected from supergravity and string theory, and are problematic in 
'tracking' quintessence models which generally have Planckian vacuum expection values 
(vevs). Couplings such as /3QF 2 /M, where F 2 is the usual Maxwell Lagrangian and M 
is the mass-scale at which we expect supergravity to fail as an effective theory cause 
variation of the fine-structure constant, and because of the large Q-vev require fine- 
tuning of the dimensionless coupling (3 of order f3 < 10~ 5 . Since we have no reason 
to expect (3 to differ significantly from order unity, this unexplained fine-tuning is 
unsettling. Carroll has argued that such dimension-five operators may be excluded 
by the existence of a discrete Z2 symmetry in the fundamental description, which acts 
on the extra dimension as <fi — > —(f), but even in this case such fine-tuning persists 
with higher-order operators of the form Q n F 2 /M n [HE]- One of our aims is to consider 
the effect of such non-renormalisable couplings in models with a low-scale of quantum 
gravity. 

There are many constraints one can consider. Despite arising from a variety of 
different physics, these bounds and the constraints on fermion couplings arising from 
stability of the classical potential share a common feature: they are sensitive to some 
power of the ratio M/M cut _ Q fj, where M is some energy scale characteristic of the process 
in question, and M cut _ g is an energy scale which controls the details of heavy physics, 
which we consider to have been integrated out in our effective description. 

In the conventional kind of four- dimensional cosmology one would usually take 
^cut-off to be of order the Planck mass Mp = G~ 2 ~ 10 19 GeV, although there are other 
natural candidates: the GUT scale M GUX around 10 15 GeV; the string scale, possibly a 
few orders of magnitude less than the Planck mass; or, more speculatively, the SUSY 
scale, which may be as low as M SUSY ~ 100 GeV. In recent years, inspired by ideas 
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from the strongly coupled limit of the heterotic string EDI ED ESI EH] an alternative 
scenario for cosmology has become popular, in which the various gauge and matter fields 
which comprise our universe are affixed to a hypersurface in a larger, five-dimensional 
bulk spacetime. This spacetime is generically a patch of Schwarzschild-anti de Sitter 
(SAdS) space. In these models, the fundamental scale Mp of quantum gravity might 
be much lower, perhaps only of order a TeV (10 _16 Mp = 10 12 eV) or so, in which case 
one would expect to obtain very different constraints on quintessence. (However one 
should note that in these models, one cannot really be dealing with the heterotic string, 
since in such theories the string scale — which amounts to the quantum gravity scale — is 
fixed to roughly coincide with the four- dimensional Planck scale and there is not much 
freedom to move it. On the other hand, theories such as Type I string theory can 
acceptably accommodate a low string scale.) Because the troubling bounds and fine- 
tunings outlined above depend sensitively on the details of the cutoff scale, one should 
carefully recalculate the constraints they impose in models with quantum gravity at low 
energies, but this is not sufficient. There are other effects which one should also take 
into account, arising from new physics associated with the branes. Most notably, for 
example, these models contain a tower of Kaluza-Klein modes in addition to a massless 
four-dimensional graviton. These Kaluza-Klein modes can be considered to arise from 
gravity in the bulk. These modes are typically massive, with masses m > 3H/2 where 
H is the brane Hubble parameter EU ^] • The presence of these modes introduces 
processes, absent in four dimensions, where quintessence can interact with gravity off 
the brane, or with the Kaluza-Klein hierarchy. 

In this paper, we apply all these ideas to constrain quintessence couplings and 
energy scales in TeV scale Planck mass models. Higher dimensional and brane-world 
models offer interesting new insights into quintessence cosmology [3H E3 UH1 1231 EE] > but 
before adopting these models wholesale it is important to consider potential constraints 
and compare them with the corresponding constraints on 4d quintessence |TH | I24 [ IT2]. 

The plan of this paper is as follows. In the next section we discuss the issue of non- 
renormalisable couplings. Then we briefly review the bounds on quintessence couplings 
in four dimensions, paying particular attention to bounds on couplings to fermions. 
In Section 0] we calculate probability amplitudes for some representative gravitational 
processes, where bulk gravitons mediate quintessence couplings to fermions. We also 
estimate the lowest-order contribution of virtual graviton exchange to the vacuum 
polarization H*(p) of the quintessence field. In Section El we state our conclusions. 
In an appendix, we give a brief derivation of the gravitational propagator in the brane 
world, using the Fadeev-Popov technique. This has appeared in the literature before 
PU] but we present this alternative derivation for simplicity and to make our account 
self-contained. 

Throughout we work in units where h — 1, but the gravitational coupling in D 
dimensions is k 2 d = 8ir/Mp~ 2 , where Mp is the D-dimensional Planck mass. We use eV 
as units of dimensionful quantities everywhere. 
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2. Non-renormalisable couplings between quintessence and the standard 
model 

One of the original motivations for TeV-scale quantum gravity was its ability to obviate 
the need for low-energy supersymmetry by removing the huge hierarchy between the 
weak and Planck scales. Such a low scale for quantum gravity can be achieved rather 
elegantly in models of large extra dimensions jU E] or in brane-world models with a 
non-compact extra dimension. 

Viewed from this perspective, the four-dimensional theory on the brane is simply an 
effective theory arising from the dimensional reduction of a more fundamental, higher- 
dimensional theory by integrating out physics at scales above which the extra dimensions 
become visible. As such, we must expect our effective four- dimensional theory to 
be burdened with a potentially infinite number of non-renormalisable interactions, 
suppressed by powers of the mass scale at which the effective theory breaks down or at 
which new physics enters the problem. In this case it is natural to assume the cutoff scale 
^cut-off to be the Planck scale of the higher- dimensional theory, M cut _ off ~ M P ~ TeV 
provided that this is lower than the effective 4d Planck scale. Above the scale M cu t-ofr 
the theory is no longer well-approximated by a four-dimensional theory. 

In the discussion that follows we will leave the cutoff mass scale M cut _ ofT arbitrary. 
To minimise clutter, we write this scale as M. Our conclusions are strongest in models 
where M ~ TeV and weaken as M increases towards the four dimensional value 
M P ~ 10 16 TeV. 

Let us begin by considering non-renormalisable Lagrangian operators of the form 
(see eg. [12]) 

Pn-£~£* (1) 

where £4 is any dimension- four standard- model operator such as F 2 or G 2 . We are 
seeking constraints on the dimensionless coupling (3, and in particular its dependence 
on the cut-off scale, as discussed above. The case £4 = F 2 leads to cosmic variation of 
the fine-structure constant a [El EE] which, assuming slow variation of Q, gives 

Aa^-n/W 1 " 1 ^ (2) 

where Q = Q(0); the symbol AQ abbreviates the field interval AQ(z) = Q(z) — Q(0); 
and z denotes redshift. The Webb et al results suggest a variation of a at the level 
Aa ^ = (— 0.543±0.116) x 10 -5 , ie. evidence for variation of a at the 4.7<r level jl^j. For 
n > 1 the ratio Q/M clearly plays a key role. We have argued that in TeV-scale higher 
dimensional models M may be close to the TeV scale. The estimate of the vev (Q) is 
model dependent; however, we can obtain constraints on standard tracker quintessence 
models where the dark energy equation of state tracks that of the dominant energy 
component of the Universe until a low redshift z ~ 1. In such models the quintessence 
field has been rolling since very early times and it is natural that it has a large vev 
today, whether the fundamental theory be four- or higher dimensional. 
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We can get a rough lower-bound on the vev of Q today by the following argument. 
To ensure that the Universe accelerates, the field Q must satisfy the standard slow-roll 
conditions familiar from early universe inflation. That is, we demand e,?) C 1 where 
e = Mp(Vq/V) 2 and rj = M^,(Vqq/V). At this point it is appropriate to comment on 
the mass scale appearing in e and rj and whether we really should be using the four- 
dimensional Planck mass Mp, or the fundamental scale M To answer this one 
can examine the effective four- dimensional equation of motion for Q. Slow-roll requires 
that the potential term Vq be sub-dominant with respect to the friction term coming 
from the Hubble expansion. At least at low energies, the effective Friedmann equation 
giving the Hubble expansion typically contains Mp, not M, and hence the effective 
four-dimensional scale Mp is the appropriate scale. 

To quantify the constraints coming from requiring e,rj <C 1, consider a standard 
tracking potential V(Q) oc Q^ 1 . In this case the field evolves as Q oc t 2 ^ 2+ ^ and the 
slow-roll parameter 77 is 



with a similar constraint arising from 77 <C 1. Clearly one can make the vev as small as 
one likes by fine-tuning the exponent 7 to be sufficiently small but then one loses the 
attractiveness of the model, since one is simply converging to a cosmo logical constant. 
Even then the fine-tuning required is significant. To have 77 < 1 with Q < TeV requires 
7 < 10~ 16 . This is quite undesirable. 

With this restriction on the form of the potential in mind, let us return to the 
dimensionless couplings (3. If we use the standard result that tracking models typically 
require vevs of order the 4d Planck scale, Q ~ Mp, then this implies that the ratio 
Q/M can be as large as Mp/TeV ~ 10 16 . Matching the Webb et al data for n = 1 is 
possible with either a fine-tuning of j3 ~ 10 -5 (the result of 12J) or is compatible with 
(3 ~ 1 by requiring AQ ~ 10~ 5 M = 10 MeV. Such a slow variation of the field is not 
consistent with the assumption of a rolling quintessence. However it is consistent with 
the Albrecht-Skordis [2] model where Q becomes trapped at the minimum of a potential 
well. Interestingly, a more detailed analysis shows that the union of constraints on a 
at various redshifts favour very little variation of Q at low redshifts z < 2 even when 



For n > 1 the situation is much worse, since the large ratio (Q/M) appears. 
The required fine-tuning on the dimensionless couplings j3 n , or variation AQ, become 
enormous. For n = 2 we have fyAQ ~ 10~ 9 eV, or if we conservatively set AQ ~ M 
we find p 2 ~ 1CT 21 . Such a coupling is quite inexplicably small. As n increases 
the fine-tuning on the dimensionless couplings (3 n rapidly increases. Of course this 
simply underlines a more fundamental point: since (Q/M) > 1 we have absolutely no 
control over the effective potential of the quintessence field. The potential should be 
computed from the higher- dimensional theory. This is a standard argument against 
chaotic inflation which typically requires super-Planckian initial conditions to obtain 
sufficient inflation and the correct amplitude of density perturbations. 




Ml 



(3) 



M = M P |35| . 
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Of course, one may simply argue that the non-perturbative potential for Q is 
completely well-behaved and only gives rise to small couplings to the standard model 
fields. This is reminiscent of the runaway dilaton model of quintessence ^Hj in which a 
massless dilaton runs to infinity where it decouples from all matter, as in the proposal 
of Damour & Polyakov [Hj. If the quintessence field is a radion, representing the 
distance between two branes, driven apart, for example, by a repulsive Casimir force, 
then we expect the Q/M — > oo limit to be trivial: the potential should vanish. This 
occurs because the Q field dynamics effectively vanish in this limit: it is equivalent to 
the degenerate case where the branes approach each other. Couplings to gauge fields 
localised on the brane world could also reasonably be expected to vanish in that limit. 
However, to appreciate this one requires the full five-dimensional picture, and in other 
scenarios the corrections may not be so harmless. 



3. Quintessence couplings in four dimensions 

In the above analysis, we obtained restrictions on a given set of coupling constants and 
shape parameters for some popular, rather generic potentials. In doing so we assumed 
that the form of the potential could simply be given as an Ansatz, so strictly speaking we 
were dealing with renormalized quantities, and the potential was the quantum effective 
potential. A more subtle question is to ask how a given tree-level potential may be 
modified when quantum effects are taken into account. This involves the study of loop 
corrections to the quintessence potential and the couplings of Q to other fields. 

For scalar quintessence and fermions, this was first done by Doran & Jackel and 
by Horvat [21], who considered couplings to neutrinos. We review their arguments as 
applied in four dimensions, and explain how this generalizes to the brane world. Many 
of the bounds described in Ref. |2Ij do not depend on the scale of gravity and are 
not strongly modified in the brane world, so we focus on the gravitational couplings 
described in Ref. [TH]. In particular, we are interested in the coupling of quintessence 
to fermions, for which the strongest bounds apply. 

In the quintessence sector we work with the Euclidean action 

S = J d 4 x V=9 (^daQ^Q + V(Q) + + m(Q))i/?j (4) 

where Q is the quintessence field, V(Q) is its classical potential, ip is a Dirac fermion, 
and m(Q) is a possibly Q-dependent fermion mass. The leading correction to the 
quintessence potential in the fermionic sector comes from the diagram of Fig. ^ 

Let m be given by a large field independent mass m plus some correction c 
generated by couplings to other fields. Then the condition that the classical potential 
dominates becomes (15] 

V L 1 AW 

T = 4^^r~ « ( 5 ) 

where we have discarded the c-independent piece m\ which does not affect dynamics, 
and we are assuming that c -C m . One can estimate V [TS] by supposing that Q 
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Figure 1. Lowest order fermion loop contributing to effective quintessence potential. 
The dashed lines represent quintessence particles, and solid lines are fermions. 




Figure 2. Low order gravitationally mediated diagrams contributing to the fermion- 
quintessence coupling. Fermions are solid lines; gravitons are wavy lines. Quintessence 
is shown as dashes, and since the graviton couples to the entire quintessence potential, 
this in principle involves arbitrary powers of Q. 



currently dominates the energy density of the universe, so that V must be comparable 
to p cr it = 8.1 x 10 -11 h 2 eV 4 . Setting A at around the GUT scale A = 10 -3 M P , and 
taking the field-independent fermion mass to be around the supersymmetry breaking 
scale, perhaps of order a TeV, or itlq = 10~ 16 Mp, that gives 

c < l(T n eV. (6) 

This calculation only depends on the details of quantum field theory in the four- 
dimensional world, so it is valid on the world volume of a brane universe provided that 
we take the effective cut-off A to be sized appropriately. Since the bound on c scales 
with A -2 , this means that a reduced cut-off will weaken any constraint. For example, in 
a model where A should be 0(TeV) ~ 10 12 eV, one finds c <C 10~ 44 eV. This weakening 
is a mixed blessing. It is harder to rule out any given quintessence model, but it may 
make the construction of a viable phenomenological model easier. 

4. Gravitationally mediated couplings in the brane world 

This bound Eq. ©, in brane or four- dimensional form, is rather stringent and could 
conceivably be violated by gravitational couplings. The low order diagrams showing 
the gravitational coupling of quintessence to fermions are shown in Fig. El In four 
dimensions, both of these diagrams involve the classical quintessence potential, so in 
fact the bound Eq. © does not apply. This happens because one can absorb the 
corrections into a renormalization of V(Q) [15J. We will compute the brane world 
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diagrams equivalent to Fig. |21 

The principal result we shall require is the propagator for gravitons in the bulk. 
As a first approximation, we work in the original Randall-Sundrum scenario where the 
branes are exactly flat. The metric is 

d^^-d^ + d^ + d, 2 ). (7) 

The propagator has been derived by Giddings, Katz & Randall 20J, and is given by 



pTsmn ^r{m ^n)s ^rs^mn (10) 



Arsrani ^ rsmn f -ik-(x-y)^ K 2 (kR) 

lon-brane (R P J k K X {kR) 1 ' 

where K v is the Macdonald or Basset function 

K u {z)= l -^ l H^\iz) ) (9) 

the brane is fixed at z = R (for RS branes, R = and p rsmn satisfies 

1 

d- 2 

where d = 8\ is the trace of the R 3 Kronecker delta. This is 3 in the present case but 
changes if one sends R 3 to R". We give an alternative derivation of this result in an 
appendix. 

It is convenient to use the Feynman rules in configuration space, rather than 
momentum space. This is because the brane matter theory has support only on the brane 
and its couplings naturally include a term 5(z—R) which is most easily accommodated in 
the configuration space formulation. The presence of the brane breaks bulk translational 
isometries in the transverse direction, and there is no conserved Noether charge to play 
the role of a conserved momentum in the z direction. For this reason, loop diagrams 
will still involve an integration over four- momenta on slices z = constant, and not full 
five-momenta in the bulk. 

The on-brane gravitational propagator is as given in Eq. ([A. 18)1 . and the fermion and 
scalar propagators are as usual. We introduce a matter theory on the brane described 
by the analogue of the four-dimensional quintessence-fermion system 



Sbrane = / d^W" det h 
Jz=R 



^d a Qd a Q + V(Q) + + m)^ 



where, h a b is the pull-back of the five dimensional metric g a b to the brane, Q is the 
quintessence field and ip is a four-dimensional (not five-dimensional) Dirac fermion. We 
are ignoring any gravitational coupling to ip via the spin connexion. In the Randall- 
Sundum case, h a b is just four-dimensional Minkowski space plus the tensor perturbation 
etj evaluated at z = R, so 

, — - — - _ tre tre 2 (tre) 2 
V^teth = l + - J- + 8^ + "' ( } 

The vertices for this theory are shown in Fig. El (See also, eg., Ref. [22] •) 

We can now proceed to evaluate the diagrams in Fig. |21 with the brane world 
graviton propagator. 
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Figure 3. Some vertices in the quantum field theory corresponding to Eq. flllfl- where 
Dirac indices are suppressed. Fermions are solid lines. Gravitons are indicated by 
wavy lines, and quintessence by dashes. The vertices are written in configuration 
space and taken to occur at a spacetime point x. Indices ij, rs on graviton lines 
refer to the R 3 index structure. Fermions entering the diagram at a point x with 
momentum p carry a coefficient function (27r) _3 / 2 u(p)e ip ' a; , where • denotes the flat, 
Euclidean inner product on the brane. Fermions leaving the diagram from a point x 
with momentum p carry {2ir)3/ 2 u(p)e~ lp ' x . Here, the constant functions u, u are Dirac 
spinors with indices suppressed: our conventions for spinors and 7-matrices match 
Weinberg ■ Quintessence particles entering or leaving the diagram with momentum 
p carry (2n)~ 3 / 2 e ±lp ' x . To find amplitudes, one integrates over the coordinates x±, 
X2, ... of all interaction points. These integrals should include the appropriate volume 
measure y — det g, which is unity on the brane. In addition, scalar products are taken in 
the metric h a b with is the pull-back of the spacetime metric g a b to the slice z = constant 
over which they are evaluated; this reduces to the flat, Euclidan scalar product on the 
brane. 



4-1. Loop diagram 

Consider the first diagram in Fig. |21 The amplitude for this process is 

• (13) 



himn 2 n( Pl )n(p 2 ) r 1 1 K 2 (kR) K 2 {k'R) 

l = TT^M y Wd) ^9/2 5 O^P) / E 77 



k'=p\— P2-k 



2 PR? K ^ cu (2tt) 9 / 2 J * kk' K l {kR)K 1 {k'R) 

where J2p = Pi ~ P2 — Pout and d 4 &E is the four-dimensional Euclidean volume measure. 

This diagram is difficult to evaluate for finite momenta p\ andp2, so we shall work in 
an approximation where all external three-momenta vanish: that is, Pi = P2 = Pout = 0- 
This approximation matches Ref. [15 j . In this case k' 2 = k 2 , and the loop integral 
becomes somewhat more tractable. In addition, in our conventions, the product of 
spinor coefficient functions u(p')u(p) evaluates to 2 when summed over spins at zero 
momentum. Temporarily ignoring the various numerical pre-factors, one has to evaluate 
the integral 



2 

Ki{kR) j 



A "*f$#§Y. (14) 
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We have explicitly written in an upper cut-off at Euclidean momenta k ~ A and a lower 
cut-off at k ~ ix. This extremely simple regularization has the advantage that it is easy 
to apply in the present context. 

The Macdonald functions K v have asymptotics governed by 

T(v) i ' z\~ v rw 

K v {z) -> ~y~ ( 2 ) as z ^ 0; K v {z) -> y —e~ z as z -> +oo. (15) 

In the infra-red, aside from numerical factors, the ratio K2(kR) / K\(kR) , behaves as a 
function of like k^ 1 . Combining this behaviour with the factor of k~ l already present 
in the propagator Eq. (|A.18|) . one can see that any infra-red divergence ought to be 
the same as in four dimensions. However, the large- z asymptotics of K v (z) changes 
the divergent behaviour in the ultra-violet. To make an estimate of Eq. f|14)) . we write 
Itx = i/i + Ii/r anc ^ approximate the integrand using its asymptotic form in both 
regions (after changing variable to z — kR), 



' d H^yj ~vl R - + Vh ^^--ln^ + -A^. 



(16) 



We have discarded a term of order 0(1/ R 2 ), which should be a good approximation 
provided A 2 3> 1/R 2 . For an extra dimension of order 1 mm, R" 1 ~ 1.97 x 1CT 4 eV, so 
this is abundantly satisfied. Eq. (fTH|) lets us pick out the leading order divergence in the 
ultra-violet and infra-red by making a small- or large-A; approximation in the integrand, 
as appropriate. 

One sees that the ultra-violet divergence, which is logarithmic in the four- 
dimensional case [IS], is modified to a considerably worse quadratic divergence. It is 
natural to interpret this modification as due to interactions with the Kaluza-Klein tower. 
Despite this, the induced coupling remains proportional to the classical quintessence 
potential V(Q), so this correction term does not destroy properties of the classical 
dynamics. This is entirely analogous to the situation in four dimensions. 



4-2. Triangle diagram 

The amplitude for the second ('triangle') diagram of Fig. |2]is 

m 2 K 4 V(Q cl ) r M(p2)H(^-ff)+mMPi) 1 K 2 (kR) 1 K 2 (k'R) 

T t 2 R 2 (2tt) 9 / 2 d[ ^ P) J E ( pi _A;) 2 + m 2 k K^kR) k> K x (k> 'R)' [ ' 

which is to be evaluated at k! = pi — p 2 — k. Using the Dirac equation, which says 
—if l u(pi) = mu(pi), the numerator can be rewritten u(p 2 )[i^ + 2m]u(pi). At zero 
external momentum, p\ — p 2 ~ and when summed over spins ^(0)^(0) = 2. One 
can also approximate ik a u(0)j a u(0) = ikoU^^u^O) = k , with a further factor of 
two arising from a sum over spins. Therefore, the numerator is just 2ko + 4m, before 
Euclidean continuation. At this point, one would typically complete the square in the 
denominator and drop terms which are odd in k a , because the integral ought to be 
rotationally invariant. However this procedure is inconvenient in the present case, since 
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Figure 4. Gravitational contribution to the quintessence self-energy. The incoming 
and outgoing particles depicted by dashes are quintessence particles. The wavy line 
circulating in the loop is the graviton. 



one wishes to keep the argument of the Macdonald functions simple. Keeping track of 
factors of i gives 

, 1R .mVTjQd) v {** (K 2 (z)V [* sin 2 g(l-cosfl) 

where we have changed variable to z = kR. This can be approximated using the same 
technique applied above for the loop diagram. We find, 

K 2 {z)\ 2 r ,„ sin 2 #(l-cos#) w , A + VA 2 + 4m 2 
" " ' 1 1 d6 — ^— -f- In 




ifi(z)/' Jo z 2 + Am 2 R 2 cos 2 6 2 A 



7t n (l + VI + 4m 2 ) 

-I In— A / fig) 

+ m 2 i? 2 /i + V/i^ + 4m 2 ' 1 ' 

where A is a reference energy scale of order m. This exhibits a logarithmic divergence in 
both the infra-red and the ultra-violet, but remains proportional to V(Q C \) and so can 
be absorbed into a redefinition of the potential. Just like the loop-diagram calculated 
above, it does not destroy the classical potential. 

This result is rather general. Since these amplitudes couple only to the quintessence 
potential through a vertex factor, which does not change when moving four dimensional 
cosmology to the brane world, the result is the same, even though the character of the 
divergences has changed. 

4-3. Gravitational contribution to quintessence mass 

As a final application of the propagator Eq. (jA.18|) . we suppose that the quintessence 
particle has some bare mass uiq and calculate the shift produced by the contribution of 
the graviton loop in Fig. 0] This vertices in the diagram take the form (— Vjrrvqdij/A and 
come from the mass term rnqQ 2 /2 in the potential and the single-graviton coupling tr e/2 
to the quintessence field. This diagram is a contribution to the self-energy i(27r) 4 n*(p) 
of the quintessence, at momentum p. That gives, with propagators for the external 
quintessence particles stripped off, 

4vrim^ 2 rAR K 2 (z) r ,„ ^ 2 sm 2 i 



W 4 n» = j fR d; j^r J o ao z2 + im2 Rlcos2e (20) 
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where we have aligned the polar axis with p, mq is the quintessence mass, we 
are assuming that the particle is on shell, and we have made the standard change 
of variable z = kR. With our conventions, the renormalized propagator becomes 
A'(p) oc (p 2 + m 2 — n*(p)) _1 , so that a negative contribution to II* gives a positive 
5m 2 . This gives the estimate 



„ m%K 2 1/2 u R , 

5ml In + -A . 21 

5 4tt 2 tR 2 \ m Q R 2m Q + ^ 2 + Am 2 Q 2 1 

That is, the divergence is linear in the ultra-violet and logarithmic at infra-red. Setting 
mq to be currently of order the Hubble rate, or mq ~ 2.1 x 10~ 33 h eV, and the infra-red 
cutoff fi to the same, we find 

5m 2 5 ~ 2.1 x KT 141 /i 3 eV 2 (5 dimensions). (22) 

This is very small, and implies that interactions with new gravitational physics 
associated with the brane world, such as the Kaluza-Klein hierarchy and graviton 
transmission through the bulk, do not seriously affect quintessence: its major problems 
remain its couplings to and interactions with normal matter. On the other hand, the 
estimate Eq. ()22|) is several orders of magnitude smaller than a comparable estimate for 
a four-dimensional cosmology with Planck scale Mp = 10 19 GeV: 

1 

8tt 

In this case, there is no infra-red divergence so the magnitude of the effect is controlled 
by the ultra-violet region. The balance between Eq. ()22|) and Eq. (J2*3*j) is controlled by 
the ratio 

5m 2 M 2 f 2 

= 4a^, (24) 
5m 2 4 Mf A K J 

where M 4 , M 5 are the four- and five- dimensional Planck scales, respectively; A is 
the four- dimensional ultra-violet cutoff; i is the five-dimensional AdS scale; and 
a = ln(l + -y/5) is a constant coming from the five- dimensional infra-red cutoff. In the 
brane world, the interpretation of this difference involves interactions with the Kaluza- 
Klein tower, whereas in the five-dimensional picture one interprets the change as a result 
of the modified Planck scale. 



5ml = _ K 2 m ^ A _ 3 Q x kt 127 ^ 3 eV 2 . (23) 



5. Conclusions 



We have studied the constraints that arise on TeV-scale quintessence models from 
a variety of sources. Non-renormalisable operators in four dimensions are typically 
important implying that the quintessence potential needs to be computed from a higher- 
dimensional framework. This follows from the fundamental mismatch between the scale 
M ~ TeV which determines the scale at which non-renormalisable operators become 
important and the vacuum expectation value, Q, of the quintessence field which is 
typically of order the 4d Planck mass in tracker models. Perturbation theory in Q/M 
fails spectacularly. 
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In contrast, the gravitational coupling of quintessence to fermionic matter in 
cosmologies of the Randall-Sundrum brane world type does not yield significant 
constraints. This is easy to understand since one expects the couplings of quintessence 
to ordinary matter to be severely constrained and sensitive to the value of the effective 
ultra-violet cut-off A UV - The brane world significantly reduces the value of this cut-off, 
and so one would expect quite radically different constraints on quintessence. 

We find that one-loop effects introduce quantum corrections in the effective 
potential just proportional to the classical potential V and therefore can just be absorbed 
into a renormalization of V. This is exactly the same as in the four-dimensional world 
and occurs for the same reason: the vertices in the diagram generate factors of V, not 
the propagator, and since this is the only quantity which changes when one moves to 
the brane world the type and character of the divergences one encounters changes, but 
the couplings remain the same. 

We have also computed the lowest-order contribution from graviton loops to the 
vacuum polarization of quintessence. In this case one must make a numerical estimate, 
and we find that the brane universe typically induces a mass shift Sm 2 very much smaller 
than in four dimensions. This shift is cut-off dependent, and scales with the ratio of the 
four- and five- dimensional Planck scales, the AdS curvature scale, and the inverse of the 
ultra-violet cutoff. From the point of view of an observer on the brane, we interpret this 
as the result of interactions with the Kaluza-Klein hierarchy and with bulk gravitons. 
The magnitude of this effect would render it undetectable and in practice the dominant 
contributions to 5m 2 would come from matter fields on the brane. 
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Appendix A. The gravitational propagator in the brane world 

In this appendix we very briefly derive the graviton propagator in the Randall-Sundrum 
scenario. Since this calculation has already appeared elsewhere j2Hj we omit details 
where they coincide. The authors of Ref. [20] deduced the propagator by solving for 
the gravitational Green's function whereas we employ the Fadeev-Popov procedure, 
but naturally the final answers shall agree. The principal result of this appendix is 
Eq. (IA.18J) for the on-brane propagator, which was used in the main text for the diagram 
calculations in Section EJ 

We adopt the conventional line element [HI Ej, for a given maximally symmetric 
three-metric 7^, 

ds 2 = -n 2 (t, y) At 2 + a 2 (t, y)^ dx i dx j + dy 2 , (A.l) 

where y is a Gaussian normal coordinate transverse to the brane. This metric is taken 
to be a solution of the five-dimensional Einstein equations with cosmological constant A 
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but vanishing bulk energy-momentum tensor. The brane is considered to be imbedded 
at y = 0, and there is a Z 2 symmetry which acts on the Gaussian normal coordinate as 
y i— > —y. There is typically a coordinate horizon where the Gaussian normal coordinates 
used in Eq. flA.l|) break down, and we write the location of this horizon as y = y^. 

Gravitational disturbances take the form of small perturbations h ab to the metric: 
ds 2 = (g ab + h ab ) dx a dx b . In a general D-dimensional spacetime, h ab will transform as 
a representation of the isometry group SO(l, D — 1). This describes the full degrees of 
freedom of the graviton. Alternatively, one could decompose h ab into its representations 
under the brane isometry group, which consists of a tensor (in the dx l dx J sector of 
the metric) and supplementary vector and scalar pieces (respectively, for vectors, in the 
dtdx\ dydx 1 sectors and for scalars in dtdy, dt 2 and dy 2 sectors) which must be added 
in to complete the full degrees of freedom of the graviton. In this paper we will deal 
only with the case of a flat, Minkowski brane which possesses a larger isometry group 
and allows us to re-absorb the vector and scalar pieces into the tensor perturbation. For 
this reason, we only calculate the tensor propagator in this section. 

This piece of the perturbation is written and takes the metric form 

dx i = -n 2 (t, y) dt 2 + a 2 (t, y)(<S y + e^) dx l dx j + dy 2 . (A.2) 
After integrating by parts and discarding surface terms, the action is 



e rs (A.3) 



S = --^ / d 5 x na 3 e^ \25 i{j 5 r]s {U - A) - ^6 l[r d 3] d s 
We have introduced an operator □ describing t and y derivatives, 

□ -— — — + — U- -) — - ( 3- + — 1 — (A 4) 

n 2 dt 2 dy 2 n 2 \ a n) dt \ a n J dy 

and A is the Sij Laplacian. The brane isometry group R 3 now appears as an invariance 
of this action: Eq. ()A.4jl is invariant under R 3 transformations (spatial rotations and 
translations). 

One passes to the quantum theory by defining correlation functions of 
the field using the functional integral prescription, (eij(x) • ■ ■ e mn (y)) = 
J[de rs ] eij(x) ■ ■ -e mn (y) expiS[e rs }. To render this integral (formally) finite, one divides 
out by the volume of the gauge group. This is the Fadeev-Popov procedure |X'6| . 
amounting to the inclusion of an extra gauge-fixing contribution in the action, 

S^S + ^r [d 5 xna 3 ^-(de b -ad b e) 2 (A.5) 

where £ and a are arbitrary numbers. We have suppressed irrelevant indices by writing 
de b = d a e ab and e = tre a b. Notice that in these formulas we are considering e ab to 
be a full five-dimensional tensor which is zero on t or y indices. Where derivatives are 
contracted with e ab this makes no difference (d a e ab = d l ei b ), but where two derivatives 
become contracted with themselves as in d a d a one must include contributions from the 
t and y sectors. 
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After integrating by parts one obtains the gauge-fixed action, 

A" 



S = -— / d 5 xna 3 e ij < 
Ik 1 J 



□ 



a 2 



1 fl_^ 

'^Oj r Oj s + OijO rs I ^ 



->-(r^) + >-(H)r (A - 6) 

In order to simplify this result, it is convenient to set £ = 1, a = 1/2 in which case both 
terms involving a disappear and one is left with the reduced action 

S = -- / d 5 x na 3 e lj D ijrs e rs (A.7) 



2 

where Dy> s is the operator 



D « M = 771 2( M*)i - ) ( □ - • ( A -8) 



These choices for a and £ coincide with the four- dimensional case. The propagator 
A rsmn (x, y) satisfies 

D tjrs A rsmn (x, y) = ~i5\x - y)6$S% (A.9) 

where 6 5 is the covariant 5-function. 

In a quite general homogeneous brane world, the metric functions n and a are not 
equal and both depend on t and y, so one has the three Killing vectors d/dx % which 
generate translations along the spacelike coordinate axes, but no other translational 
Killing vectors. For this reason, it is sensible to try and diagonalize A rsmn as a Fourier 
transform in the x l , but one will not be able to deal with the t and y dependence in the 
same way. Therefore, 

A rsmn = _. p rsmn J _^P_ e -ip-(«»-v) G (p ; XQj y Q] x ^ y ^ ( A .10) 

where p IS the combination 

d 2 

and d = 8\ is the trace of the R 3 Kronecker delta. We are now adopting a convention 
of writing the coordinates of any point x in spacetime as (a; , a?, a; 5 ). Substituting 
Eq. (jA.lOjl into Eq. (|A.9J) shows that G must obey the equation 



□ 



^ }G = 2k 2 5(x° -y°)5(x 5 -y 5 ). (A.12) 
a 2 J 



Appendix A.l. The Randall-Sundrum propagator 

At this point, one can make no further progress without specifying some explicit form 
for a and n. The simplest choice is to take the brane to be empty of matter, except for 
some intrinsic tension which is tuned to give a Minkowski brane. The line element 
is 

ds 2 = e~ 2%l (-dt 2 + S i:j dx l dxJ) + dy 2 . (A.13) 



Radiative constraints on brane quintessence 



17 



In this special case, the functions a and n do turn out to be equal, and all quantities 
are independent of the cosmic time t, so one recovers d/dt as a Killing symmetry. 
This is a great convenience, because one can now write G as a Fourier integral 

in x° — y°, leaving only an ordinary differential equation 
for the x 5 , y 5 dependence of the Fourier transform G. Changing to a conformal bulk 
coordinate z defined by dx 5 = adx z , this ordinary differential equation turns out to be 
just the Bessel equation, 

BG = -2^j^5(x z -y z ), (A.14) 

where B is the Bessel operator, 

B = + + (/3 2 - . (A.15) 

d(x z f x z dx z V ( x ") J 

In these coordinates, the location of the brane is z = but we will often take its 
location to be arbitrary and write z = R instead. When making numerical estimates, 
we restore IR = 1. We have introduced a new quantity j3 defined by f3 2 = u 2 — p 2 and 
define a four-vector k = (u,p) with k 2 = —f3 2 in our signature diag (—1, 1, 1, 1). 

From this point, out derivation coincides with the earlier derivation of Giddings, 
Katz & Randall 20 , so we omit further details and merely state the result. The general 
solution for G is a combination of Bessel functions, and by integrating Eq. (JA.14J) over 
a small neighbourhood of x z = y z one obtains a continuity condition on G and a step 
condition on dG/dx z . Demanding that the normal derivative of G vanish at the brane, 
and that positive frequency waves be purely ingoing in the far-field, together with the 
junction conditions at the brane allows one to solve uniquely for G. The result for the 
entire propagator is 

d 3 p did 



A rsmn = p rsmn J _^_^_ e - ip .( x -y) +kl ;( x - y °) yy(x*, y Z ] R), (A.16) 

where W(x z ,y z ; R) satisfies 



W{x z ^R) = -i(-\ {Ji{PR)bP(Pz<) ~ H?>{pR)W*>)) (A.17) 

\y j 1 ti x \pR) 

in which z = R is the location of the brane and respectively min{x z ,y z }, 

max{x z ,y z }. is the Hankel function of order v of the first kind. This propagator 
agrees with the expressions (2.14)-(2.15) given in Ref. j20j once notational differences 
have been taken into account. 

There is an obvious simplification of Eq. (jA.16|) in the special case that both 
endpoints x, y are taken on the brane. One finds, 

lK 2 {kR) 

£R P J (2vr) 4 kK^kR) 
where K v is the Macdonald or Basset function 

K u {z) = ~7rf +1 ifW(iz) (A. 19) 



a rsmn I _ lK „rsmn f a K -jfc.fo-y) 1 A 2 \Ktt) , . 

^ lon-brano — e r> P I ^_N4 e i. r/ /i„n\ {A.LO) 
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and is entirely real. 

In the text, we make use of Wick rotation to Euclidean signature. Since the graviton 
propagator has acquired a much more complicated space-time dependence than the 
propagators of Minkowski space fields, there may be more obstructions to this procedure 
than can be dealt with by moving the l/k poles of Eq. ifPSjl off- axis. In particular, 
if the Macdonald functions K v have singularities anywhere on C, then a Wick rotation 
could not be justified. To see that this is not so, it is convenient to make use of the 
following integral representation of the Macdonald function 

KJz) = <[f - , . fe-'r 1 / 2 ( 1 - M "~ 1/2 df (A.20) 
V ; V 2zT(v + 1/2) Jo V 2zJ K ' 

From this it is clear that K v has a singularity at z = for any v but is otherwise analytic 
everywhere. Moreover it is not zero except at z = oo. There is an essential singularity 
at oo stemming from the exponential, but since the contour remains fixed at z = and 
z = oo, this does not interfere with the analytic continuation. 
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